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The treatment of the electron— nucleus interaction based on the Mott dif- 
ferential cross section was extended to account for effects due to screened 
Coulomb potentials, finite sizes and finite rest masses of nuclei for electrons 
above 200 keV and up to ultra high energies. This treatment allows one to 
determine both the total and differential cross sections, thus, subsequently to 
calculate the resulting nuclear and non-ionizing stopping powers. Above a few 
hundreds of MeV, neglecting the effect due to finite rest masses of recoil nu- 
clei the stopping power and NIEL result to be largely underestimated. While, 
above a few tens of MeV, the finite size of the nuclear target prevents a further 
large increase of stopping powers which approach almost constant values. 



1. Introduction 

Nuclei and electrons populate the heliosphere. Most of the nuclei are galac- 
tic cosmic rays (GCR), while electrons can additionally be originated by 
the Sun and Jupiter's magnetosphere, which is a major source of relativistic 
electrons in the heliosphere (e.g., see Ref.^'^ and references therein). Pro- 
tons and electrons are also major constituents of the Earth's radiation 
belts. These particles can interact with materials and onboard electro- 
nics in spacecrafts, inducing displacements of atomic nuclei, thus inflicting 
permanent damages. As the particle energy increases, for instance above 
« 20 MeV for protons and « 130 MeV/nucleon for a-particles (e.g., see 
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Section 4.2.1.4 and Figure 4.26 at page 418 of Ref.''), the dominant me- 
chanism for displacement damage is determined by hadronic interactions; 
for electrons and low-energy nuclei the elastic Coulomb scattering is the 
relevant physical process to induce permanent damage. 

The non-ionizing energy-loss (NIEL) is the energy lost from particles 
traversing a unit length of a medium through physical processes resulting 
in permanent atomic displacements. The displacement damage is mostly 
responsible for the degradation of semiconductor devices - like those using 
silicon - where, for instance, depleted layers are required for normal ope- 
ration conditions (e.g. see Ref.^). The nuclear stopping power and NIEL 
deposition - due to elastic Coulomb scatterings - from protons, light- and 
heavy-ions traversing an absorber were previously dealt"'''*' with (sec also 
Sections 1.6, 1.6.1, 2.1.4-2.1.4.2, 4.2.1.6 of Ref.^). In the present work, the 
nuclear stopping power and NIEL deposition due to elastic Coulomb scat- 
terings of electrons are treated up to ultra relativistic energies. 

The developed model (i.e., see Sects. 2-2.4) for screened Coulomb elastic 
scattering up to relativistic energies is included into Geant4 distribution'' 
and is available with Geant4 version 9.5 (December 2011). In Sects. 3, 4, 
the nuclear and non-ionizing stopping powers for electrons in materials are 
treated, while a final discussion is found in Sect. 5. 

2. Scattering Cross Section of Electrons on Nuclei 

The scattering of electrons by unscreened atomic nuclei was treated by 
Mott* (see also Sections 4-4.5 in Chapter IX of Ref.^) extending a method 
of Wentzel^° (see also Born^^) and including effects related to the spin of 
electrons^ . Wentzel's method was dealing with incident and scattered waves 
on point-like nuclei. The differential cross section (DCS) - the so-called 
Mott differential cross section (MDCS) - was expressed by Mott^ as two 
conditionally convergent infinite series in terms of Legendre expansions. In 
Mott-Wentzel treatment, the scattering occurs on a field of force generating 
a radially dependent Coulomb - unscreened (screened) in Afott^ (Wentzcl-'^") 
- potential. It has to be remarked that Mott's treatment of collisions of fast 
electrons with atoms (e.g., see Chapter XVI of Ref.^) involves the knowledge 
of the wave function of the atom, thus, in most cases the computation 
of cross sections depends on the application of numerical methods (see a 
further discussion in Sect. 2.2). Furthermore, the MDCS was derived in the 
laboratory reference system for infinitely heavy nuclei initially at rest with 
negligible spin effects and must be numerically evaluated for any specific 
nuclear target. Effects related to the recoil and finite rest mass of the target 
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nucleus (M) were neglected. Thus, in this framework the total energy of 
electrons has to be smaller or much smaller than Mc^ . 

As discussed by Idoeta and Legarda^^ (e.g., see also Refs.^^'^*), Mott 
provided an "exact" differential cross section because no Born approxi- 
mation^ of any order is employed in its derivation. Various authors have 
approximated the MDCS for special situations, usually expressing their re- 
sults in terms of ratios, 1Z, of the so-obtained approximated differential cross 
sections with respect to that one for a Rutherford scattering (RDCS) - the 
so-called Rutherford's formula, see Section 1.6.1 of Ref.^ - for an incoming 
particle with z =1 given by: 

da^^^ fZe'V 1 f Ze^\'^ 1 



dO. \pl3cj (l-cos6»)^ VPpcJ sin^(6l/2) 

1 



^2mc2/32-yy sin^(6l/2)' 

where m is the electron rest mass, Z is the atomic number of the target 
nucleus, (3 ~ v/c with v the electron velocity and c the speed of light; 7 is 
the corresponding Lorentz factor; p and are the momentum and scattering 
angle of the electron, respectively; finally, since the interaction is isotropic 
with respect to the azimuthal angle, it is worth noting that dQ can be given 
as 

dQ = 2'n:sin0d0. (2) 
The MDCS is usually expressed as: 

dci ~ dn ' ^ ' 

where 7^"^°** (as above mentioned) is the ratio between the MDCS and 
RDCS. In particular, Bartlett-Watson-*^^ determined cross sections for nu- 
clei with atomic number Z = 80 and energies from 0.024 up to 1.7 MeV (see 
also Ref.^^). McKinley and Feshbach^^ expanded Mott's series in terms of 
power series in aZ (with a the fine-structure constant) and {aZ)/(3] these 
expansions, which give results accurate to 1% up to atomic numbers Z k AQ 
(e.g., see discussions in Refs.^*'^^), were further simplified to obtain an ap- 
proximate analytical formula with that accuracy for aZ < 0.2. Feshbach^*^ 
tabulated values of the differential cross section as a function of scattering 



''In quantum mechanical potential scattering, the scattered wave may be obtained from 
the so-called Born expansion. The Born approximation is the first term of the Born 
expansion (see, for instance, references indicated in Section 1.6.1 Ref.^). 
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angle for nuclei with atomic number up to 80 and electrons with kinetic 
energies larger than 4 MeV. Curr^* reported values of the differential cross 
section as a function of scattering angle accurate at 1% for {aZ)/ j5 < 0.6; 
while Doggett and Spencer^^ tabulated the MDCS for energies from 10 
down to 0.05 MeV. Recently, Idoeta and Legarda^^ provided a further se- 
ries transformations and made a systematic comparison with those from 
McKinley and Feshbach^'' , Curr^^ , Doggett and Spencer^^ . For electrons 
with kinetic energies from several keV up to 900 MeV and target nuclei with 
1 ^ ^ 90, Lijian, Quing and Zhengming^^ provided a practical interpo- 
lated expression [Eq. (16)] for 7^"^°** with an average error less than 1%; in 
the present treatment, that expression - discussed in Sect. 2.1 - is the one 
assumed for 7^^^'!°** in Eq. (3) hereafter. 

The analytical expression derived by McKinley and Feshbach^^ - men- 
tioned above - for the ratio with respect to Rutherford's formula [Equa- 
tion (7) of Ref.^^] is given by: 

^McF = 1 _ ^2 gjjj2 (^/2) + Z a/37r sin(e/2) [1 - sin(^/2)] (4) 
with the corresponding differential cross section (McFDCS) 

where da^^^/dV, is from Eq. (1). It has to be remarked that for positrons, 
the ratio TZ^o^ becomes 

nf^f = 1 - /?2 sin2(6'/2) - Za/37rsin(6'/2) [1 - sin(0/2)] (6) 

(e.g., see Equation (6) of Ref.^^ ). Furthermore, for Mc^ much larger than 
the total energy of incoming electron energies the distinction between la- 
boratory (i.e., the system in which the target particle is initially at rest) 
and center-of-mass (CoM) systems disappears (e.g., see discussion in Sec- 
tion 1.6.1 of Ref.'^). Furthermore, in the CoM of the reaction the energy 
transferred from an electron to a nucleus initially at rest in the laboratory 
system (i.e., its recoil kinetic energy T) is related to the maximum energy 
transferable Tmax as 

T = T^^ sm^{9'l2) (7) 

[e.g., sec Equations (1.27, 1.95) at page 11 and 31, respectively, of Ref.^], 
where 0' is the scattering angle in the CoM system. From Eqs. (2, 7) one 
obtains 

dT = d^l'. (8) 

477 ^ ' 
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Since 9 is ^ 9' for M c? much larger than the electron energy, one finds that 
Eq. (7) can be approximated as 

T~T^3,sin2(0/2), 



and 



sin2(^/2) = 



47r 



(9) 
(10) 

(11) 



Using Eqs. (4, 10, 11), Eqs. (1, 5) can be respectively rewritten as: 



" dT 



4tt dT 



rp2 ' 



(12) 



rp2 



T 



+ Za/3TT 



T 



T 



PI3c 



rp2 



1-/3- — (^+Za7r) + Za/3m 



(13) 



Ze^YnT, 



pj3c 



max ^McF 



rp2 



with 



1 -1^7^ (/3 + Za-n:) + Za/37ri 



(14) 



[e.g., see Equation (11.4) of Ref.,^'' see also Ref.^^ and references 
therein]. Similarly, for positrons one finds 



da. 



McF 



Zf\ ttT, 



rji2 



1 - /9- — (^ - Za-n) - Za^n. 

max 



T 



[e.g., sec Rcfs.^^'^^ and references therein]. Finally, in a similar way the 
MDCS [Eq. (3)] is 



dT 



dT 



\ '^'^max rpMott (rp\ 

-30 T2 



(15) 



December 7, 2011 1:55 



WSPC - Proceedings Trim Size: 9in x 6in Consolandi2011 



with 7e'^°"(T) from Eq. (18). 

2.1. Interpolated Expression for T^^ott 

As mentioned in Sect. 2, Curr^^ derived TZ^^°*-^ as a function the atomic 
number Z of the target nucleus and velocity /3c of the incoming electron at 
several scattering angles from ^ = 30 up to 180 . Recently, Lijian, Quing 
and Zhcngming^^ provided a practical interpolated expression [Eq. (16)] 
which is a function of both 9 and /3 for electron energies from several keV 
up to 900 MeV, i.e., 

4 

^Mott =J2aj{Z,p){l - cos^)j/^ (16) 
j=o 

where 

6 
k=l 

and pc = 0.7181287c is the mean velocity of electrons within the above 
mentioned energy range. The coefficients 6k,j(-^) are listed in Table 1 
of Ref.22 for 1 ^ Z 90. 

At 10, 100 and 1000 MeV for Li, Si, Fe and Pb, values of 7^'^°" were 
calculated using both Curr-*^^ and Lijian, Quing and Zhengming^^ meth- 
ods and found to be in a very good agreement. It has to be remarked 
that with respect to the values of Tl^'''^ obtained from Eq. (4) at 100 MeV 
one finds an average variation of about 0.2%, 3.2% and 8.8% for Li, Si 
and Fc nuclei, respectively. However, the stopping power determined using 
Eq. (52) (i.e., with 7^'^°") differs by less than 0.5% with that calculated 
using Eq. (53) (i.e., with Tl^"^). Tl^"^^ obtained from Eq. (16) at 100 MeV 
is shown in Fig. 1 for Li, Si, Fc and Pb miclci as a function of the scattering 
angle. Furthermore, it has to be pointed out that the energy dependence 
of T^^'^o** from Eq. (16) was studied and observed to be negligible above 
« lOMoV [as expected from Eq. (17)]. 

Finally, from Eqs.(7, 16) [e.g., see also Equation (1.93) at page 31 
of Ref.'^], one finds that 7^^!°'-'' can be expressed in terms of the transferred 
energy T as 

4 / \j/2 

n^^'\T) = Y,a,{Z,P)(^] . (18) 

j=0 V-imax/ 
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2.2. Screened Coulomb Potentials 

As already mentioned in Sect. 2, a complete treatment of electron intera- 
ctions with atoms (e.g., see Chapter XVI of Ref.^) involves the knowledge of 
the wave function of the target atom and, thus - as remarked by Fernandez- 
Vera, Mayol and Salvat^^ -, a relevant amount of numerical work when the 
kinetic energies of electrons exceed a few hundreds of keV. 

The simple scattering model due to Wentzel^^ - with a single expo- 
nential screening function [e.g., see Equation (2.71) at page 95 of Ref.^ , 
Equation (21) in Ref.^^ and Ref.-'^^] - was repeatedly employed in treat- 
ing single and multiple Coulomb scattering with screened potentials 
(e.g, see Ref.^^ - and references therein - for a survey of such a topic 
and also Refs.^'^'^^"^*). Neglecting effects like those related to spin and 
finite size of nuclei, for proton and nucleus interactions with nuclei it was 
shown that the resulting elastic differential cross section of a projectile 
with bare nuclear-charge ez on a target with bare nuclear-charge eZ dif- 
fers from the Rutherford differential cross section (RDCS) by an additional 
term - the so-called screening parameter - which prevents the divergence 
of the cross section when the angle 9 of scattered particles approaches 
0° [e.g., see Refs.^'^'^^"^^ (see also references therein) and Section 1.6.1 
of Ref.'^]. It has to be remarked that the RDCS for z = 1 particles can also 
be employed to describe the scattering of non-relativistic electrons with 
unscreened nuclei (e.g, see Refs.^'^^ and references therein). As derived by 
Moliere^^ for the single Coulomb scattering using a Thomas-Fermi poten- 
tial, for ^; = 1 particles the screening parameter Ag^u [e.g., see Equation (21) 
of Bethe^''] is expressed as 



where a, c and h are the fine-structure constant, speed of light and reduced 
Planck constant, respectively; p {(3c) is the momentum (velocity) of the 
incoming particle undergoing the scattering onto a target supposed to be 
initially at rest - i.e., in the laboratory system -; axF is the screening length 
suggested by Thomas-Fermi (e.g., sec Rcfs.^^''^^) 




(19) 



axF 



ZV3 



(20) 



with 
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the Bohr radius, m the electron rest mass and 

Ctf = 2 ( ^ ) - 0.88534 

a constant introduced in the Thonias-Fcrmi model [e.g., see Equa- 
tions (2.73. 2,82) - at page 95 and 99, respectively - of Ref.^ and 
Ref.^ , see also references therein]. The modified Rutherford's formula 
[da^^ {9) / (Kl], i.e., the differential cross section - obtained from the 
Wentzel-Moliere treatment of the single scattering on screened nuclear po- 
tentials - is given by [e.g., see Equation (2.84) of Ref.^ , Section 2.3 in Ref.^^ 
and Ref.^ (see also references therein)]: 



da^^'{0) _ f zZe 



dfl \2ppcJ [A,,M + sin2(^/2)]' 

_ (i(7^"^ sin'^(^/2) 
~ [^,M + sin2(^/2)]' 



(21) 



dn 



r{0). (22) 



with 



_ sin^(g/2) 
^^'^ A,M+sin^(^/2)- 

^(0) - the so-called screening factor - depends on the scattering angle 6 
and screening parameter Ag^u- As discussed in Sect. 2.4, in the DCS the 
term Ag^u cannot be neglected [Eq. (22)] for scattering angles (6) within a 
forward (with respect to the electron direction) angular region narrowing 
with increasing energy from several degrees (for high-Z material) at 200 keV 
down to less than or much less than a mrad above 200 MeV. 

An approximated description of elastic interactions of electrons with 
screened Coulomb fields of nuclei can be obtained factorizing the 
MDCS, i.e., involving Rutherford's formula [da^^^/d^] for particles with 
z = 1, the screening factor ^{9) and the ratio 7^^^°** between RDCS and 
MDCS: 



dfl dQ 



S^{9) 7^'''^°" (24) 



[e.g., see Equation (1) of Ref.-^^ , Equation (A34) at page 208 of Ref.-^^ , see 
also Ref.^^ and citations from these references]. Thus, the corresponding 
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Fig. 1. -RMott obtained from Eq. (16) at 100 MeV for Li, Si, Fe and Pb nuclei as a 
function of scattering angle. 



screened differential cross section derived using the analytical expression 
from McKinley and Feshbach^^ can be approximated with 

It has to be remarked - as derived by Zeitler and Olsen^^ - that spin and 
screening effects can be separately treated for small scattering angles; while 
at large angles (i.e., at large momentum transfer), the factorization is well 
suited under the condition that 

2Z4/V-i-«l 

(e.g., see Refs.^^''^^). Zeitler and Olsen'^^ suggested that for electron ener- 
gies above 200 keV the overlap of spin and screening effects is small for 
all elements and for all energies; for lower energies the overlapping of the 
spin and screening effects may be appreciable for heavy elements and large 
angles. 
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2.3. Finite Nuclear Size 

As suggested by Fernandez- Vera, Mayol and Salvat^^ , above lOMeV the 
effect of the finite nuclear size has to be taken into account in the treat- 
ment of the electron-nucleus elastic scattering. With increasing energies, 
deviations from a point-like behavior (see, for instance. Figure 4 of Ref.,"'^^ 
j^g£ 32,33 g^j^^ references therein) were observed at large angles where the 
screening factor [Eq. (23)] is ~ 1. 

The ratio between the actual measured and that expected from the 
point-like differential cross section (e.g., the MDCS) expresses the square 
of the nuclear form factor {\F\) which, in turn, depends on the momentum 
transfer q, i.e., that acquired by the target initially at rest: 

VT(T + 2Mc^) 
c 

with T from Eq. (7) or, for Mc^ larger or much larger than the elec- 
tron energy, from its approximate expression Eq. (9) [e.g., see Equa- 
tions (31, 57, 58) of Ref.33 , Section 3.1.2 of Ref.^ , Refs}^'^^'^^'^% 

The factorized differential cross section for clastic interactions of elec- 
trons with screened Coulomb fields of nuclei [Eq. (24)] accounting for the 
effects due to the finite nuclear size is given by: 



dfl dfl 



Fiq)\ 



2 



J Rut 

^ ^ ^'^^^ ^"""^ '^^^^l' ^^^^ 

[e.g., see Equation (18) of Ref.^** , Ref.-^^ and also references therein]. Thus, 
using the analytical expression derived by McKinley and Feshbach^^ 
[Eq. (4)] one obtains the corresponding screened differential cross section 
[Eq. (25)] accounting for the finite nuclear size effects, i.e., 



dn dn 

S''{e)n'^'=''\F{q)\'' (28) 



^^Rut ^McF I 7?/„m2 



dn 

X {l-/32sin2(6'/2)-t-Za/37rsin(6l/2) [1 - sin(6'/2)]} . (29) 
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In terms of kinetic energy, one can respectively rewrite Eqs. (27, 28) as 

J Mott (rp\ j^Rut 

^^'^ ^ ' t{T) n^°'\T) \F{q)\' (30) 



rf<f (T) da^-\T) ^McF 



dT dT 

dT - dT ^""'^(^^ l^^^^l' 

with da^'^ydT from Eq. (12), TZ^°*\T) from Eq. (18), n^''^{T) from 
Eq. (14) and, using Eqs. (7, 9, 23), 

T 

The nuclear form factor accounts for the spatial distribution of charge 
density probed in the electron-nucleus scattering [e.g., see Equation (58) 
of Rcf.^^ , Section 3.1.2 of Rcf.^ , Rcfs.^'''^^'^^'^'* and references therein]. For 
instance, among those spherically symmetric treated in literature, one finds 
that for i) an exponential charge distribution (i^exp) [e.g., see Equation (6) 
of Ref.^* , Equation (93) at page 252 of Ref.^'^ and references therein], ii) 
a Gaussian charge distribution (i^gau) [e.g., see Equation (6) of Rcf.^*^ and 
references therein] and iii) an uniform-uniform folded charge distribution 
over spheres with different radii {F^) [e.g., see Equation (22) of Ref.^^ , 
Ref.^^ and references therein]. For instance, the form factor i^exp is 



^'exp(g) 



1 



12 \ h ) 



(32) 



where is the nuclear radius [e.g., see Equation (6) of Ref.^*]. To a first 
approximation, rn can be parameterized by 

r„ = 1.27^"-2^ fm (33) 

with A the atomic weight [e.g., see Equation (7) of Ref.^®]. Equation (33) 
provides values of rn in agreement up to heavy nuclei (like Pb and U) 
with those available, for instance, in Table 1 of Ref.^^ . The nuclear form 
factor is 1 for (7 = and rapidly decreases with increasing q [e.g., see 
Eq. (32), Equation (6) of Ref.^^ and Equation (22) of Ref." for Fexp, Fg^u 
and Fu, respectively]. Furthermore, from inspection of Eqs. (7, 9, 26) small 
q are those corresponding to scattering angles within the forward (with 
respect to the electron direction) angular region which, in turn, narrows 
with increasing electron energy. For instance, in lithium the square values 
{\F{q)f) of these form factors are in agreement within 1% up to 6' < 124.1° 
(2.4°) at 20MeV (IGeV); in silicon up to 9' < 138.4° (2.4°) at 20MeV 
(IGeV); in iron up to 9' < 108.0° (2.1°) at 20MeV (IGeV). However, as 
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Silicon 



1 1 I 1 I 1 1 I I 

10° 10^ 10' 10^ 10" 10= 10^ 

Kinetic energy [MeV] 

Fig. 2. As a function of the kinetic energy of electrons from 200 keV up to 1 TeV, ratios 
of nuclear stopping power of electrons in silicon calculated neglecting i) nuclear size effects 
(i.e., for |Fexp|'^ = 1) (dashed curve) and ii) effects due to the finite rest mass of the 
target nucleus (dashed and dotted curve) [i.e., in Eq. (53) replacing da^^'^p q^j^{T) / dT 
with da^'^^{T)/dT from Eq. (31)] both divided by that one obtained using Eq. (53). 

discussed in Sect. 2.4, these upper angles are larger or much larger with 
respect to those required to obtain 99% of the total cross section. Thus, 
the usage of any of the above mentioned nuclear form factors - e.g., Fcxp as 
in the present treatment - is expected to be appropriate in the treatment 
of the transport of electrons in matter, when single scattering mechanisms 
are relevant, for instance in dealing with the nuclear stopping power and 
non-ionization energy-loss deposition. 

2.4. Finite Rest Mass of Target Nucleus 

The DCS treated in Sects. 2-2.3 is based on the extension of the MDCS 
to include effects due to interactions on screened Coulomb potentials of 
nuclei and their finite size. However, in the treatment, the electron energies 
were assumed to be small (or much smaller) with respect to that (Mc^) 
corresponding to the rest mass (M) of target nuclei. 

The Rutherford scattering on screened Coulomb fields - i.e., under the 
action of a central force - by massive charged particles at energies larger or 
much larger than Mc^ was treated by Boschini et al.^'^ in the CoM system 
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(e.g., sec also Sections 1.6, 1.6.1, 2.1.4.2 of Rcf.'' and references therein). It 
was shown that the differential cross section [da^^^{9')/d^' with 9' the 
scattering angle in the CoM system] is that one derived for describing the 
interaction on a fixed scattering center of a particle with i) momentum 
p[ equal to the momentum of the incoming particle (i.e., the electron in 
the present treatment) in the CoM system and ii) rest mass equal to the 
relativistic reduced mass /Xrei [e.g., see Equations (1.80, 1.81) at page 28 
of Ref.^] . /irei is given by 

A*rel = Tl— (34) 



mMc 



(35) 



where p is the momentum of the incoming particle (the electron in the 
present treatment) in the laboratory system: m is the rest mass of the 
incoming particle (i.e., the electron rest mass); finally, M\^i is the invariant 
mass - e.g.. Section 1.3.2 of Rcf.'' - of the two-particle system. Thus, the 
velocity of the interacting particle is 



0,c = c 



\ 



v. 



(36) 



[e.g., see Equation (1.82) at page 29 of Ref.^]. For an incoming particle with 

2 = 1, da^^\d')ld^ is given by 



with 



dQ! 



'^P'rPrCj [A,+sin2(0'/2)] 



2 ' 



2p; OTF 



1.13 + 3.76 X 



aZ 



the screening factor [e.g., see Equations (2.87, 2.6 
of Ref.^] . Equation (37) can be rewritten as 



da^^'{e') _ da^^'^'ie') ^2 



with 



dn' dQ' 



2p'J',cj sin*(0V2) 



(37) 



(38) 



at page 103 



(39) 



(40) 
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the corresponding RDCS for the reaction in the CoM system [e.g., see 
Equation (1.79) at page 28 of Ref.^] and 

^°°^^^^-A + sin^(^'/2) ^ ^ 

the screening factor. Using, Eqs. (7, 8), one can respectively rewrite 
Eqs. (40, 41, 39, 37) as 



(42) 

i?CoM(T) = ^— - (43) 



5m(T) = — 

^ max^ 



dT dT 



dcoM{T) (44) 



da'^^'iT) (Ze'^X^ T, 



dT ""UP^cJ (T^^^A+Tf ^^^^ 

[e.g., see Equation (2.90) at page 103 of Ref.^ or Equation (13) of Ref.^]. 

As already mentioned (Sect. 2.2), the screening parameter Ag prevents 
the DCS to diverge - see last term in Eq. (37) -, i.e., for 9' of the order of 
or smaller than 



6*^^ = arcsin (2^/A^ 



effects due to screening of the nuclear Coulomb field have to be accounted 
for. 6'^^ rapidly decreases with increasing the kinetic energies of elec- 
trons. For instance, in iron 6'^^ is « 1.7° (0.03 rad) at 200 keV and w 0.004° 
(7.0 X 10-2 mrad) at 200 MeV; in silicon, it is w 1.3° (0.022 rad) at 200 keV 
and w 0.003° (5.5 x 10"^ mrad) at 200 McV; while, in lithium, it is w 0.75° 
(13mrad) at 200 keV and w 0.002° (3.3x10-^ mrad) at 200 MeV. Therefore, 
in Eq. (39) the term Ag (i.e., the screening parameter [Eq. (38)]) cannot be 
neglected for scattering angles within a forward angular region narrowing 
with increasing energies from a few degrees (for \ow-Z material) at about 
200 keV down to less than or much less than a mrad above 200 MeV. It is 
worthwhile to remark that in silicon, for instance, 6' can be approximated 
with 6 up to a few hundred McV. 

To account for the finite rest mass of target nuclei, the factorized MDCS 
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[Eq. (27)] has to be re-expressed in the CoM system as: 



^^'^ SLuiO') «(^') \F{q)\\ (46) 



dQ' 

where F{q) is the nuclear form factor (Sect. 2.3) with q the momentum 
transfer to the recoil nucleus [Eq. (26)]; finally, as discussed in Sect. 2.1, 
•^Mott exhibits almost no dependence on electron energy above w 10 MeV, 
thus, since at low energies i2 9' and /3 ±2 7e^°M (^') is obtained replacing 
9 and with 9' and respectively, in Eq. (16). 

Using the analytical expression derived by McKinley and Feshbach^^ , 
one finds that the corresponding screened differential cross section account- 
ing for the finite nuclear size effects [Eqs. (28, 29)] can be re-expressed as 

dn' — dn/ — ^comk^) 7^coM(^') \Fm (47) 

with 

^Com(^') = {1-/3,' sin2(^'/2) + Za^>sin(^72) [l-sin(^72)]} . (48) 

It has to be remarked that scattered electrons are mostly found in the 
forward or very forward direction. For instance, using Eq. (48) one can 
derive that in lithium f« 99% of electrons are scattered with 6' < 0.27° 
(0.007°) at 20 MeV (IGeV); in silicon with 6' < 0.46° (0.009°) at 20 MeV 

(IGeV); in iron with 9' < 0.6° (0.013°) at 20 MeV (IGeV). 

In terms of kinetic energy T, from Eqs. (7, 8) one can respectively rewrite 
Eqs. (46, 47) as 



J Mott (rp\ j^Rut' 

^ = ^ dkuiT) «(T) |F(,)|' (49) 



dT dT 



tcouiT) 7^S(^) \F{q)\' (50) 



with da^^'^'/dT from Eq. (42), 5coM(r) from Eq. (43) and n^^l^{T) re- 
placing P with in Eq. (14), i.e.. 



J- max V -i T 



(51) 



December 7, 2011 1:55 



WSPC - Proceedings Trim Size: 9in x 6in Consolandi2011 



16 



5.0x10 




= 2.0x10 



10 10 
Kinetic Energy [MeV] 



10° 



Fig. 3. Nuclear stopping powers (in MeVcm^/g) in '^Li, ^^C, ^*Si and ^®Fe - calculated 
from Eq. (53) - and divided by the density of the material as a function of the kinetic 
energy of electrons from 200 keV up to 1 TeV. 

Finally, as discussed in Sect. 2.1, 7?.'^°**(T) exhibits almost no dependence 

on electron energy above « 10 MeV, thus, since at low energies d' and 
/3 ^ 'R^^°m{T) is obtained replacing /3 with /J^ in Eq. (18). 



3. Nuclear Stopping Power of Electrons 

Using Eq. (49), the nuclear stopping power - in MeV cm" 
electron-nucleus interaction can be obtainc^d as 

Mott „T, 



dE 
dx 



dT 



TdT 



of Coulomb 



(52) 



nucl 

with ua the number of nuclei (atoms) per unit of volume [e.g., see Equa- 
tion (1.71) of Ref ^] and, finally, the negative sign indicates that energy 
is lost by electrons (thus, achieved by recoil targets). Using the analytical 
approximation derived by McKinley and Feshbach,^'' i.e., Eq. (50), for the 
nuclear stopping power one finds 

McF .Tm 



-f-Y 



= nA 

Jo 



^<tcoM(r) 

dT 



TdT. 



(53) 



As already mentioned in Sect. 2.4, the large momentum transfers - cor- 
responding to large scattering angles - are disfavored by effects due to 
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10° 10^ 10^ 10' 10" 10' 10° 
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Fig. 4. For = 21 eV, ratios of NIELs of electrons in silicon calculated as a function of 

the kinetic energy of electrons from 220 keV up to 1 TeV neglecting i) nuclear size effects 
(i.e., for |FexpP = 1) (dashed curve) and ii) effects due to the finite rest mass of the 
target nucleus (dashed and dotted curve) [i.e., in Eq. (55) replacing da^'^p ^^j^{T) / dT 
with da^'^^{T)/dT from Eq. (31)] both divided by that one obtained using Eq. (55). 



the finite nuclear size accounted for by means of the nuclear form fac- 
tor (Sect. 2. 3). For instance, in Fig. 2 the ratios of nuclear stopping powers 
of electrons in silicon arc shown as a function of the kinetic energies of 
electrons from 200 keV up to 1 TeV. These ratios are the nuclear stopping 
powers calculated neglecting i) nuclear size effects (i.e., for |Fexp|^ = 1) and 
ii) effects due to the finite rest mass of the target nucleus [i.e., in Eq. (53) 
replacing da^^lcouiT) / dT with da'^§{T)/dT from Eq. (31)] both divided 
by that one obtained using Eq. (53). Above a few tens of MeV, a larger 
stopping power is found assuming |i^exp|^ = 1 and, in addition, above a few 
hundreds of MeV the stopping power largely decreases when effects due to 
the finite nuclear rest mass are not accounted for. 

In Fig. 3 , the nuclear stopping powers in ''Li, ^^C, ^^Si and ^^Fe are 
shown as a function of the kinetic energy of electrons from 200 kcV up to 
ITeV. These nuclear stopping powers in MeVcm^/g arc calculated from 
Eq. (53) and divided by the density of the medium. The flattening of the 
high energy behavior of the curves is mostly due to the nuclear form factor 
which prevents the stopping power to increase with increasing T„iax- As 
expected, the stopping power are slightly (not exceeding a few percent) 
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Fig. 5. For = 21 eV, NIEL (in McV cm'^ g^""^) in silicon calculated using Eq. (55) as a 
function of the kinetic energy from 220 keV up to 1 TeV (dashed and dotted curve) ; NIEL 
values from Messenger et al.*^ (o) and Jun et al.*^ (•) calculated in the laboratory system 
without accounting for the effects due to the screened Coulomb potential, finite size and 
rest mass of recoil silicon; the dotted curve is obtained replacing da^'f. ^^y^{T) / dT with 
da^''^{T)/dT [Eq. (13)] in Eq. (55). 



varied at large energies replacing Fexp with i^gau or (Sect. 2.3). However, 
a further study is needed to determine a most suited parametrization of the 
nuclear form factor^^"^^ particularly for high-Z materials; for instance, in 
lead the stopping power results to be depressed at energies of about (20- 
40) MeV, while in medium and light nuclei this occurs at energies of the 
order or above 100 MeV. 



4. Non-Ionizing Energy-Loss of Electrons 

A relevant process - which causes permanent damage to the silicon bulk 
structure - is the so-called displacement damage (e.g., see Chapter 4 of Ref.,^ 
Refs.*'^'^^ and references therein). Displacement damage may be inflicted 
when a primary knocked-on atom (PKA) is generated. The interstitial atom 
and relative vacancy are termed Frenkel-pair (FP). In turn, the displaced 
atom may have sufficient energy to migrate inside the lattice and - by fur- 
ther collisions - can displace other atoms as in a collision cascade. This 
displacement process modifies the bulk characteristics of the device and 
causes its degradation. The total number of FPs can be estimated calcu- 
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lating the energy density deposited from displacement processes. In turn, 
this energy density is related to the non-ionizing energy loss (NIEL), i.e., 
the energy per unit path lost by the incident particle due to displacement 
processes. 

In case of Coulomb scattering of electrons on nuclei, the non-ionizing 
energy-loss can be calculated using (as discussed in Sect. 2-3) the MDCRS 
or its approximate expression McFDCS [e.g., Eqs. (49, 50), respectively], 
once the screened Coulomb fields, finite sizes and rest masses of nuclei are 
accounted for, i.e., in MeV/cm 

-f^) =n. / TLiT) "---y;M(T) 



or 



ydxJ^^McF Jt, dT 



(55) 



[e.g., see Equation (4.113) at page 402 and, in addition. Sections 4.2.1- 
4.2.1.2 of Ref.^], where T is the kinetic energy transferred to the tar- 
get nucleus, L{T) is the fraction of T deposited by means of displace- 
ment processes. The Lindhard partition function, L{T), can be approxi- 
mated using the so-called Norgett-Robintson-Torrens expression [e.g., see 
Refs.^^'^'o and/or Equations (4.121, 4.123) at pages 404 and 405, respecti- 
vely, of Ref.^ (see also references therein)]. T^e = T L{T) is the so-called 
damage energy, i.e., the energy deposited by a recoil nucleus with kinetic 
energy T via displacement damages inside the medium. In Eqs. (54, 55) the 
integral is computed from the minimum energy - the so-called threshold 
energy for displacement, i.e., that energy necessary to displace the atom 
from its lattice position - up to the maximum energy Tmax that can be 
transferred during a single collision process. For instance, Td is about 21 eV 
in silicon (e.g., see Table 1 in Ref.^^ and references therein) requiring elec- 
trons with kinetic energies above w 220 keV [e.g., see Equation (4.142) at 
page 412 of Ref.^]. 

As already discussed with respect to nuclear stopping powers in Sect. 3, 
the large momentum transfers (corresponding to large scattering angles) 
are disfavored by effects due to the finite nuclear size accounted for by the 
nuclear form factor. For instance, in Fig. 4 the ratios of NIELs for electrons 
in silicon are shown as a function of the kinetic energy of electrons from 
220 keV up to 1 TeV. These ratios are the NIELs calculated neglecting i) 
nuclear size effects (i.e., for l-Foxpl^ = 1) and ii) effects due to the finite rest 
mass of the target nucleus [i.e., in Eq. (55) replacing da^'^p Q^yi{T) / dT with 
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daf^'§{T)/dT from Eq. (31)] both divided by that one (Fig. 5) obtained 
using Eq. (55). Above lOMeV, the NIEL is « 20% larger assuming 
|Fexp|^ = 1 and, in addition, above (100-200) MeV the calculated NIEL 
largely decreases when the effects of nuclear rest mass are not accounted 
for. Finally, it has to be remarked that similar results can be obtained 
neglecting the screening factor: already at energies lower that 200 keV, k, 
21 eV is much larger than Tmax^s- 

In Fig. 5, for Td = 21 cV the non-ionizing energy loss (in MeVcm^ g""'^) 
calculated using Eq. (55) in silicon is shown (dashed and dotted curve) as 
a function of the kinetic energy from 220 keV up to 1 TeV and is compared 
with that one tabulated by Messenger et al.*^ (Jun et al.^^) from « 240 keV 
up to 200 MeV (1 GeV). For the laboratory system. Messenger et al.**^ used 
the approximate MDCS found by Doggett-Spencer-^^ and Lindhard's par- 
tition function numerically obtained by Doran*^ without accounting for 
the effects due to screened Coulomb potential (i.e., = 1), finite size 
(i.e., = 1) and finite rest mass of the silicon target; while, Jun et al.'''^ 
followed the approach discussed in Ref.^-^ (see the treatment in Sect. 2.1) to 
determine an approximate expression of the MDCS and dealt the Lindhard 
partition function using the modified Norgett-Robintson-Torrens expres- 
sion found*^ by Akkerman and Barak (2006). The dotted curve is obtained 
replacing da^'^^couiT) / dT with da^''^{T)/dT [Eq. (13)] in Eq. (55): at 
lOOMeV-lGeV, the agreement between the latter calculation and values 
from Messenger et al.''^ and Jun et al.^'^ is within several percents. It has 
to be remarked (see also Fig. 4) that i) above (100-200) MeV effects due 
to screened Coulomb potentials, finite sizes and finite rest masses of nuclei 
have to be taken into account and ii) for energies between w 100 MeV and 
« 1 GeV the effects of neglecting the nuclear form factor and finite rest 
mass of nuclei almost compensate each other. 

5. Conclusions 

The treatment of electron-nucleus interactions accounting for effects due 
to screened Coulomb potentials, finite sizes and finite rest masses of nuclei 
allows one to determine both the total and differential cross sections, thus, 
to calculate the resulting nuclear and non-ionizing stopping powers from 
low (about 200 keV) up to very high energy (ITeV). 



''Jun ct al. (2009) determined that the usage of the Norgett-Robintson— Torrens expres- 
sion or, alternatively, the one modified by Akkerman and BaraJj (2006) yields similar 
NIEL values. 
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Above a few hundreds of MeV, neglecting the effects of finite rest masses 
of recoil nuclei, the stopping power and NIEL result to be largely underesti- 
mated. Above a few tens of MeV the finite size of the nuclear target prevents 
a further increase of both stopping power and NIEL, which approach almost 
constant values. The fiattcning of the high energy behavior of the nuclear 
and non-ionizing energy-losses is mostly due to the nuclear form factor 
which prevents stopping powers to increase with increasing Tmax- However, 
a further study is needed to determine a most siiited parametrization of 
the nuclear form factor able to provide a satisfactory trend in the energy 
region below about hundred MeV also for high-Z materials. 

Finally, at lOOMeV-1 GeV an agreement to within several percents was 
obtained between the present calculation with respect to the NIEL values 
from Messenger et al.^'^ and Jun et al.''"^ . 
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